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Abstract 


The present paper describes a new method of de- 
sign decomposition for structural analysis and optimiza- 
tion. For this method, the structure is divided into sub- 
structures where each substructure has its structural re- 
sponse described by a structural-response subproblem, 
and its structural sizing determined from a structural- 
sizing subproblem. The structural responses of sub- 
structures that have rigid body modes when separated 
from the remainder of the structure are further decom- 
posed into displacements that have no rigid body com- 
ponents, and a set of rigid body modes. The structural- 
response subproblems are linked together through forces 
determined within a structural-sizing coordination sub- 
problem which also determines the magnitude of any 
rigid body displacements. Structural-sizing subprob- 
lems having constraints local to the substructures are 
linked together through penalty terms that are deter- 
mined by a structural-sizing coordination subproblem. 
All the substructure structural-response subproblems 
are totally decoupled from each other, as are all the sub- 
structure structural-sizing subproblems, resulting in the 
significant potential for use of parallel solution methods 
for these subproblems. 


an optimum design problem is whether any solution ex- 
ists to the decomposed problem. If a solution exists, 
another concern is whether this solution to the decom- 
posed problem is also a solution to the original prob- 
lem. The present paper addresses these concerns, and 
develops a formal method for solving structural-design 
problems which are decomposed using substructuring. 

Specifically, derivations are presented for a decom- 
position of the minimum weight design process when 
each design constraint depends only on the design vari- 
ables and structural response of a single substructure. 
Stress and local buckling constraints and certain dis- 
placement constraints fall into this category of de- 
sign constraint. For this decomposition method, the 
structure is partitioned into substructures in a man- 
ner that allows for solution of the structural response 
and the structural optimization on a substructure- by- 
substructure basis. The present decomposition method 
also has several unique aspects that allow for efficient 
computation, much of which can be done in parallel. 
The theory of the decomposition method is discussed 
utilizing a two-substructure example and the derivation 
for an arbitrary number of substructures is presented in 
the appendix. 


Introduction 

Nonlinear mathematical programming (NLP) is ex- 
tensively used for optimal structural sizing in struc- 
tural design. The most computationally efficient NLP 
methods at present incorporate approximation con- 
cepts, such as developed in reference 1, that require 
both the structural response of the entire structure and 
the response sensitivity derivatives with respect to all 
the design variables. These approximation-based NLP 
structural-sizing methods have been very successful and 
are implemented in several commercial structural anal- 
ysis codes. 

Further improvements in the computational effi- 
ciency of NLP structural-design methods are desired, 
especially for the detail design of structures requiring a 
large number of design variables and constraints. One 
approach to the solution of large problems is to decom- 
pose a large problem into a set of smaller subproblems. 
If most of the smaller subproblems can be solved sepa- 
rately and in parallel, the time required to find the opti- 
mum design can be significantly reduced. By repeating 
this decomposition process at the subproblem level, a 
multilevel decomposition (ref. 2 and references therein) 
can be obtained. One concern in the decomposition of 


Theory 

The method of determining the structural response 
of a structure that is decomposed into substructures is 
presented first. This response solution method is then 
utilized in formulating a set of Equilibrium Program- 
ming (EP) subproblems that determine the optimal de- 
sign. The derivation which follows demonstrates the 
method using two substructures which are assumed to 
have a linear elastic response. The appendix summa- 
rizes the method for an arbitrary number of substruc- 
tures. 

Decomposition of the Structural Response 

The present structural-response decomposition 
method has some similarity to the work presented in 
reference 3 in that the present method is a substruc- 
turing approach which utilizes Lagrange multipliers to 
enforce compatible displacements at the interfaces be- 
tween substructures. Four salient features characterize 
the present method. The structural response of each 
substructure is decomposed into the rigid body mo- 
tions (referred to as “modes” herein) of the substruc- 
ture, and displacements that have strain energy but no 
rigid body motion. Secondly, an augmented stiffness 



matrix is formed for each substructure. These stiff- 
ness matrices are symmetric, and can be factored in- 
dependently and, computationally, in parallel. Thirdly, 
a structural-response coordination problem determines 
the internal forces between the substructures (i.e., the 
Lagrange multipliers), and the rigid body modes. The 
structural-response coordination problem requires the 
factored, augmented stiffness matrices of the substruc- 
tures, and results in a system of linear equations with 
the number of degrees- of- freedom equal to the num- 
ber of shared degrees-of- freedom between substructures 
plus the total number of substructure rigid body modes. 
Lastly, once the internal forces between the substruc- 
tures are determined, the displacements having strain 
energy can be determined. 

The simple wing structure finite element model in 
figure 1 is utilized to demonstrate the derivation of the 
present decomposition method. The model shown in 
the figure is decomposed into two substructures, the 
second of which has rigid body modes. Only one load- 
ing condition is considered in the derivation; however, 
the derivation is easily extended to multiple loading 
conditions. The displacements of the entire wing for 
a given loading condition are denoted by the nodal dis- 
placement vector U, and the displacements of the two 
substructures are given by the vectors Ui and U2. The 
displacements at the interface of the two substructures 
must be equal for the two substructures to be compat- 
ible. Those displacements at the interfaces of the two 
substructures that must be compatible are chosen in a 
predetermined order with signed Boolean matrices de- 
noted by Bi and B2, where the nonzero entries of Bi 
are equal to -hi and the nonzero entries of B2 are equal 
to - 1 . The constraints of compatible displacements at 
the interface between substructures is then expressed by 
B1U1 -h B2U2 = 0 . The external loading for the entire 
wing is given by the nodal force vector F, and the de- 
composition of the external loading is given by vectors 
Fi and F2 for the two substructures. The decomposi- 
tion of the externally applied nodal forces at the inter- 
face nodes is arbitrary as long as the sum of the forces 
applied to the interface nodes of each substructure is 
equal to the actual externally applied forces. 

The linear elastic structural response for the entire 
wing is determined using a minimum potential energy 
formulation. This response is the solution to the un- 
constrained minimization problem given by 

min Qu t KU - F r u) ( 1 ) 

where the first term is the strain energy of the structure, 
and the second term is the work done by the external 
forces. The necessary conditions for the unconstrained 
minimization problem represented by statement ( 1 ) 


simplify to the linear equations: 

KU = F ( 2 ) 

where K is the global stiffness matrix. In a substruc- 
tural decomposition, a stiffness matrix can be formu- 
lated for each substructure. These stiffness matrices 
are relative to the substructure nodal displacement vec- 
tors Uj, and are denoted by K^. Thus a constrained 
minimization problem that is equivalent to the problem 
given by statement ( 1 ) and uses the substructure nodal 
displacement vectors is 

A G u > ' KiUi + 

-FfUi-FjU 2 ) < 3 ) 

BiUi +B2U2 = 0 

At first glance, the constrained minimization prob- 
lem given by statement ( 3 ) appears to be separable, i.e., 
the problem can be decomposed into two distinct min- 
imization problems, with design variables Ui and U2, 
respectively, and a coordination subproblem. However, 
the minimization problem for U2 would be insoluble 
due to the rank deficiency of matrix K2. Since the final 
form of the design decomposition method of the present 
paper will be an EP problem, another argument against 
the (Ui,U 2 ) decomposition can be made based on the 
sufficient conditions for the existence of an EP solution 
(see ref. 4 ). The existence theorem for an EP solution 
requires that the solution space for the EP subproblem 

Uf’K 2 U 2 -FjU 2 ) ) where F2 here includes 

interface forces from substructure 1 , be bounded. How- 
ever U2 is not bounded as seen from the following argu- 
ment. Because substructure 2 has six rigid body modes, 
a rigid body matrix R2 can be formed that consists of 
the six rigid body modes as column vectors. This rigid 
body matrix satisfies the condition K2R2 = [ 0 ], thus 
the vector R2CX2 added to any solution to this subprob- 
lem will not change the potential energy if the vector 
0L2 satisfies the scalar equation F^R2a2 = 0 . Since 
the vector a 2 is not otherwise restricted, it (and there- 
fore U2) is unbounded, thus the conditions for solution 
existence cannot be satisfied. 

From the preceding arguments, any decomposition 
must explicitly account for the rigid body modes of 
the substructures. A decomposition that does this 
follows. Let U2 = U2 + R2CK2 with the supplementary 
condition R^U2 = 0 . The minimization problem given 
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by statement (3) then becomes 

min fiufKiUi + Ju£k 2 u 2 

(Ui, 112,0:2) \2 2 

-Ff U 1 - F^(u 2 + R2Q!2 ) J 

BlUx 4 - B 2 (U 2 + R 2 a 2 ) = 0 

rJ U 2 = 0 


which are solved for Ui and (u 2 ,/x), respectfully. The 
structural-response coordination EP subproblem that 
determines A and a 2 is obtained by substituting the 
solutions of equations (7) and (8) into the last two 
equations in necessary conditions (5) 


M 

' A ' 

_ 

B 1 Kf 1 Fi+ 

[B 2 j[0]]K 2 - 1 [F!> 7 ’] T 


.“2. 


-R-2 F 2 


where the condition K 2 R 2 = [0] has been utilized, 
and the supplementary condition has been added to 
the problem given by statement (3). The necessary 
conditions for the structural-response problem given by 
statement (4) are 

KiUx -Fi +BfA = 0 


where the matrix M is given by 


M = 

BiK^B^f j 

[b 2 j [ 0 ] ]k 2 1 [b 2 j [0] ] T 1 

— b 2 r 2 


-(B 2 R 2 ) t 

[°] 


“u 2 - 


f 2 -b^a‘ 

. M - 


0 


r^(f 2 -b£a) = 0 

B 1 U 1 4 - B 2 (u 2 + R 2 Q 2 ) = 0 


( 5 ) 


where A are the Lagrange multipliers for the first con- 
straint in statement (4), /u are the Lagrange multipliers 
for the second constraint in statement (4), and the ma- 
trix K 2 is the augmented stiffness matrix given by 


K 2 ! R 2 


k 2 


( 6 ) 


Ro j [o] 


L iV 2 


Because the matrices B* in necessary conditions (5) 
are Boolean, the Lagrange multipliers A are simply 
the interface forces between the substructures. Using 
the second equation in necessary conditions (5), the 
condition R^K 2 = [0], said the fact that the columns 
of R 2 are linearly independent, the third equation 
in necessary conditions (5) implies that the Lagrange 
multipliers fi are equal to zero. 

Two structural-response EP subproblems that axe 
implied by necessary conditions (5) are those with nec- 
essary conditions 


KiUi=Fi-BfA (7) 


and 


U 2 * 


f 2 — bJa 

. M . 


0 


( 8 ) 


The EP subproblem with necessary conditions given 
by equation (8) includes the constraint that u 2 has no 
rigid body modes. This constraint essentially bounds 
the solution space, and thus satisfies the requirement 
for existence of a solution. The existence of a unique 
solution to equation (8) is guaranteed from the linear 
independence of the columns of matrix K 2 which results 
from augmenting the stiffness matrix K 2 with the rigid 
body modes R 2 that are in its null space. Thus, the 
structural response (i.e., nodal displacements) can be 
solved using the following steps: 1) factor the matrices 
Ki and K 2 that are used in equations (7) and (8); 2) use 
these factored matrices to formulate the matrix M and 
the right hand side of equation (9); 3) solve equation (9) 
for A and a 2 ; and 4) solve equations (7) and (8) for Ui 
and u 2 , respectively. 

Structural Optimization with Constraints Local 
to the Substructures 

In this section, the structural optimization problem 
is decomposed into a set of EP subproblems that are 
used to perform the optimization, and the related struc- 
tural analysis, of each substructure. Only structural- 
sizing design variables, herein termed sizing variables, 
such as plate thickness or beam section properties, are 
considered in the present derivation. Also, the present 
derivation only considers constraints local to the sub- 
structures, that is, a constraint relation must explic- 
itly depend on the sizing variables and structural re- 
sponse of only one substructure. The two- substructure 
wing example from the previous section will be used 
to demonstrate the present method of decomposition; 
however, this method is extended to multiple substruc- 
tures in the appendix. The decomposition is derived 
from the necessary conditions of a simultaneous anal- 
ysis and design formulation (ref. 5). In the present 
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simultaneous analysis and design formulation, mini- 
mum weight is the design goal, both the structural 
displacements and the sizing variables are utilized as 
design variables, and the optimization constraints are 
the equations governing structural response (treated as 
equality constraints) as well as the usual inequality con- 
straints that ensure that the design meets the strength, 
buckling, and other design requirements. The present 
method has the following features: 1) the optimization 
of each substructure is independent and can be done in 
parallel (although the overall optimization procedure is 
still iterative); 2) the resulting design is optimal; 3) the 
objective function of each substructure optimization is 
the weight of the substructure plus a penalty term de- 
rived from an structural-sizing coordination subprob- 
iem; 4) the structural-sizing coordination subproblem 
is a matrix equation that uses the same matrix M as in 
the structural-response coordination subproblem (equa- 
tion (10)); and 5) each substructure optimization de- 
pends directly on only the sizing variables and displace- 
ments of the substructure being optimized. 


The local constraint function vectors for substruc- 
tures 1 and 2 will be denoted as gi(vi,Ui) and 
g2( v 2* U2), respectively. The vectors \\ and V2 are 
vectors of sizing variables for the two substructures. In 
general, the constraint function vector for the second 
substructure is expressed by the relation g2(v2,U2) = 
g2( v 2> U2 + R2C*2)> but, because stress is invariant with 
respect to rigid body motions, it simplifies to g2(v2, U2) 
for local strength and buckling constraints. The ap- 
proach to forming the decomposed optimization prob- 
lem is as follows: 1) express the optimal design prob- 
lem as a simultaneous analysis and design problem us- 
ing necessary conditions (5) for the previously derived 
structural-response method as additional constraints; 
2) determine the necessary condition equations for this 
simultaneous analysis and design problem; 3) reduce the 
necessary condition equations for the sizing variables of 
each substructure to a form which has a contribution 
from the local substructure and a term representing 
coupling to other substructures; and 4) treat these re- 
duced necessary condition equations as if they are from 
an EP problem, and transform the reduced necessary 
condition equations to a set of EP subproblems and a 
structural-sizing coordination subproblem. 


problem to be solved is given by 


min W 1 (vi)+W 2 (v 2 ) 

(vi,V2,Ui,U2,Q2,A,/i) 

gl(vi,Ui) < 0 
g2(v2,u 2 + R.2<*2) < 0 
KiCvOUi -Fx + BfA = 0 

K 2 (v 2 ) 

B 1 U 1 + B2(u 2 4- R2a 2 ) = 0 

fi = 0 


U2 


'f 2 -b|a 

. M 


0 


( 11 ) 


where the weights W\ and and the stiffness ma- 
trices of the two substructures are explicit functions of 
the sizing variables as shown, and the third equation in 
necessary conditions (5) is replaced with the equivalent 
condition \i — 0. 


Necessary Conditions . The Lagrangian func- 
tion for this optimization problem is 


L =R r 1 (v 1 ) + W 2 (v 2 ) +7i’gi(v 1 ,Ui) 

+ g2(v2, u 2 + R.2£*2) 

+ 1 '1 (KjCvOUx-Fj+BfA) 


+ 14 K 2 (v 2 ) 


"u 2 " 


f 2 -b^a' 

- M . 


0 


+ 6 t (B1U1 + B 2 (u 2 + 1*20:2)) + r ]2 H 


( 12 ) 


Simultaneous Analysis and Design For- 
mulation . The optimal design problem is posed using 
the simultaneous analysis and design approach where 
necessary conditions (5) that describe the structural 
response are formulated as constraints. Thus, the 


In this definition, the quantities 7 j, 6 , and 77 are La- 

grange multipliers for the constraints in statement (11). 
Taking the derivatives of the Lagrangian function de- 
fined in equation (12) with respect to the design vari- 
ables of the optimization problem in statement (11) 
gives the following equations which are part of the set 
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of necessary conditions 


equation in necessary conditions (13) 


T dgi . T . T d KiUj T 

71 ^ + 1 "aTT” 0 


91^2 T<9g2 , r 9K 2 [u 2 jo j T 

■^7 + 72 ^ 2 ~ “° 


7i 


£jo r l T 

dv 2 

1 


T ^ gl +^K 1 +tf T B 1 =o t 


9Ui 


^§g+*^ 


jpxp 


[ 0 ] 


+ 6 1 B 2 = 0 J 


(13) 


72’|^R2 + 6 T B 2 R 2 = 0 T 

du2 


,.T nT I ..T 

^1 t>i + v 2 


v\ Tk 2 


[ 0 ] 


idqxq 


rn^ 1 1 


[ 0 ] 


+ V2 = 0 T 


where / r xr is the identity matrix of dimension r, p 
is the number of degrees-of- freedom for the structural 
response of substructure 2, q is the number of rigid body 
modes for substructure 2, and the relation pi — 0 has 
been utilized. 

The sensitivity relations for equations (7) and (8) 
with respect to the sizing variables will be necessary for 
continuing the derivation. These relations are expressed 

by 


dKiUi 

<9vi 


dV 


dX 


+ Kl ^ + B ^ 


[o] 


aK 2 [u7jo T ] T 

9v 2 


+ k 2 


' d U2* 
dv 2 


r B r4A i 

2 dv 2 

~ToT" 


= [0] 


(14) 

Using equations (14), the sensitivity derivatives of the 
structural response holding the interactions between the 
substructures fixed (i.e., fixed A) are 


dUi 
dv i 


K 


-i dK l V l 

dv\ 


T 


7i 


dgi 

dUi 


+ « J Bi Kj 


r l 


(16) 


and the Lagrange multipliers v<i are obtained by com- 
bining the fourth and last equations in necessary con- 
ditions (13) 


T 

"2 = " 


j r,* 


K 


-l 


(17) 


Substituting equations (16) and (17) into the first and 
second equations of necessary conditions (13), respec- 
tively, and utilizing definitions (15), gives two of the 
reduced necessary conditions 


dW i 

ctvi 


dw 2 

dv 2 


T 

7i 


+ 


dgi dUx 


aui dvi 


+6 1 Bx 


dUi 


dvi 


+ 


T ( 9 te dg2 du2 \ 

\9v 2 9u 2 dv 2 \) 


(18) 


+6 t B 2 ^ 
dv 2 


X dv 2 


= O' 


The final step in reducing the necessary condition equa- 
tions is to find expressions for <5 and 772 * Substituting 
equations (16) and (17) into the sixth equation in nec- 
essary conditions (13), and defining the quantity f ) 2 = 
— [r£r 2 ] 1 172 which implies that — 77 J [Rj j [0] ] = 
V 2 1 [0] j Iqxq ] Kj the following expression is obtained 


111 ! 
1 

.■>-,aK 2 [uf!on T (15) 

' 6 ' 

T 

BxK^Bj + [B 2 ! [0]]K2 1 [B 2 j [0]] T ‘ 

l^L 

2 ch/2 



-(B 2 R 2 ) t 

r i r 1 T 


Although the total derivative d/x/dv 2 is identically 
equal to zero as indicated in equations (14), the re- 
stricted derivative dfi/dv 2 \\ in definitions (15) is not. 

Reduction of Necessary Conditions. Pro- 
ceeding with the reduction of the necessary conditions, 
the Lagrange multipliers are obtained from the third 




k; 


b 2 j [0] 




(19) 

Combining this equation with the fifth equation in 
necessary conditions (13) gives the following equation 
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for 6 and f } 2 


equations (18), and is given by the statement 



The matrix M in equation (20) is the same matrix 
defined in equation (10). 

Transformation to EP Subproblems. The 
reduced necessary condition equations (18) can be 
transformed into EP subproblems having these equa- 
tions as their necessary conditions. Four types of 
subproblems result: structural-response subproblems, 
structural-sizing subproblems, a structural-response co- 
ordination subproblem, and a structural-sizing coor- 
dination subproblem. Only the structural-sizing sub- 
problems are actually transformed into constrained- 
minimization problems. The first step in these trans- 
formations is to define approximate models for Ui and 
U2 that do not consider coupling of the substructures. 
To fulfill the EP necessary conditions, the approximate 
models for Ui and U2 must have the same values and 
the same sensitivity derivatives as the exact responses 
for the optimal values of the sizing variables. To solve 
the individual subproblems efficiently, only simple mod- 
els that are explicit functions of the subproblem sizing 
variables are considered. Although many approximate 
models are possible (see ref. 6), a simple first-order Tay- 
lor series is used in the present derivation. The approx- 
imate models are defined by 


Uf (vj) = — - (v! - vO + Ui 

dw 1 A 

u 2 ( v 2) = ( v 2 — V 2 ) + U2 

dv 2 A 


(21) 


where the displacements with the caret accents and the 
restricted derivatives in equations (21) are calculated 
at specific values of the sizing variables denoted by 
v*. Also for compactness, the dependence on these 
previously calculated values are not explicitly given in 
the function arguments. With these definitions, two 
EP structural-sizing subproblems can be formed that 
have reduced necessary condition equations (18) as their 
necessary conditions. The EP subproblem for the first 
substructure is developed from the first equation in 


min Wi(vi) + 6 r BiU^(vi) 

( 22 ) 

gl(vi,U^(vi)) < 0 

where the Lagrange multipliers for gi are 71. Simi- 
larly, the structural-sizing subproblem for the second 
substructure is developed from the second equation in 
reduced necessary conditions (18), and is given by the 
statement 

min W 2 (v 2 ) + 6 T E^Uo 1 ^) + Wo -7^- v 2 
V2 dv 2 A (23) 

g2(v2, u^(v 2 ) + R2a=2) < 0 

where the Lagrange multipliers for g2 are 72. Equa- 
tion (20) is the coordination subproblem that deter- 
mines the penalty coefficients 6 and r )2 = — R^R2 
by utilizing the Lagrange multipliers determined from 
solving the subproblems given by statements (22) and 
(23). The structural response is given by solving equa- 
tions (7), (8), and (9), and the restricted derivatives are 
given by definitions (15). 

Discussion 

Summarizing the salient features of the present 
method, a set of interacting EP subproblems, each as- 
sociated with a particular substructure except for two 
coordination subproblems, is iteratively solved until 
the equilibrium solution is achieved. The structural- 
response subproblem for each substructure is indepen- 
dent of the structural response of the rest of the struc- 
ture because the coupling of the structural-response 
subproblems is achieved using a set of Lagrange mul- 
tipliers which represent forces between the substruc- 
tures. A computational advantage is obtained be- 
cause the expensive factorization of matrices can be 
done in parallel for each substructure. The Lagrange 
multipliers representing the forces between substruc- 
tures, and the substructure rigid body modes are de- 
termined within a structural-response coordination sub- 
problem. The optimization of substructures that have 
constraints that are local to a substructure is accom- 
plished in structural-sizing subproblems in a manner so 
that the individual substructure optimizations are in- 
dependent of each other, and can be accomplished in 
parallel. The computational efficiency is enhanced be- 
cause the only structural-response sensitivity derivative 
information required is local to the substructure being 
optimized. The penalty terms in each structural-sizing 
subproblem that represent the coupling of the substruc- 
ture optimizations ensure that the overall design is op- 
timal. Additional computational advantages may be 
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accrued because the calculation of the coefficients for 
these penalty terms is done in a structural-sizing coor- 
dination subproblem that utilizes the same matrix uti- 
lized in the structural-response coordination subprob- 
lem. The approximate updated sensitivity derivative 
method reported in reference 7 may prove to enhance 
further the computational efficiency of the method. 

Several aspects of the decomposition method require 
further research. One aspect that needs attention is the 
translation of the theory derived in the present report 
into a viable numerical algorithm. This translation im- 
plies the need for methods to form and factor the aug- 
mented stiffness matrices of the substructures, and to 
form and factor the coupling matrix M that is nec- 
essary for calculating the substructure interface forces 
and rigid body modes, and the penalty terms of the 
structural-sizing subproblems. In addition, although 
several EP theorems may prove existence of a solution, 
in practice the substructural optimizations will require 
rational selections of move limits on the sizing vari- 
ables, and a strategy for dealing with subproblems hav- 
ing no feasible solution (a modification of the method 
described in reference 6 that includes a constraint vio- 
lation penalty for each structural-sizing subproblem is 
presently envisioned). Also, the stability of the solu- 
tion process, and the convergence to the equilibrium 
solution needs to be studied. Another area of inter- 
est is the extension of the method to constraints that 
are more general than constraints which are local to a 
substructure. The utilization of parallel computers to 
implement the method also requires investigation. 

Conclusions 

A new method for structural analysis and structural 
optimization utilizing a substructure-based decomposi- 
tion has been derived in detail for two substructures and 
extended for multiple substructures. Because the start- 
ing point for the decomposition is a simultaneous design 
and analysis formulation of the design problem, the so- 
lution to the resulting decomposed system of subprob- 
lems is optimal. The final form of the design decompo- 
sition is a set of loosely coupled subproblems that de- 
fine an Equilibrium Programming (EP) problem. Thus, 
existence of a solution to the decomposed system can 
be studied utilizing EP existence theorems. All of the 
structural-response subproblems for the substructures 
cam be solved independently of each other, as can all of 
the structural-sizing subproblems. Thus, significant po- 
tential exists for the development of a parallel solution 
method that utilizes this decomposition of the optimal 
design problem to reduce the solution time for large 
problems. 



Appendix 


Extension of Method to Multiple Substructures 

Assume that there are n substructures, and that m substructures (0 < m < n) have rigid body modes when 
separated from the rest of the structure. The n* displacement degrees-of-freedom of substructure i are denoted by 
Uj = (Un , . . . , Ui n . ) T . The substructures having rigid body modes are ordered to be the last m substructures. These 
substructures have displacements denoted by U* = u;-bR;a{ where R* is the matrix containing the rigid body modes 
of substructure t, and Kj u; = 0. The equivalencing of degrees-of-freedom at common interface nodes in the different 
substructures is represented by the set A of n 1 4-tuples defined by A = { (j, k, p y q) | j < k and Uj p = U }. In 
the definition of set A, redundancies in the equivalencing of degrees-of-freedom are omitted. For example, a degree- 
of-freedom shared by three substructures need only be equivalenced between two pairs of substructures, not all 
three pair-wise combinations. The compatibility constraint equations (i.e., the equations that enforce compatibility 
between the substructures) are defined by assigning an order to the elements of set A that corresponds to the order 
of the compatibility constraint equations. Thus, the r th compatibility constraint equation will depend on degrees- 
of-freedom in substructure j if either the first or second element of the r th 4-tuple in A is j . These compatibility 
constraint equations are expressed explicitly by defining the signed Boolean matrices Bj for j = 1, . . . , n which have 
dimensions n 1 x nj. Matrix Bj has a 1 at location rp if the r th 4-tuple in A has the components (j, •), and it 

has a -1 at location rq if the r th 4-tuple in A has the components (♦, j, •,<?). Otherwise, the entries in matrix Bj 
are zero. Thus, the compatibility constraint equations for all the substructures are symbolically represented by the 
system of n 1 equations: 

n 

£BiUi = 0 (Xl) 

i= 1 

The definition of the Boolean matrix B ; in the present section reduces to the previous definition of the Boolean 
matrix for the case of only two substructures. 


Decomposition of the Structural Response 

Utilizing the previously described definitions, the minimization problem for the structural response given by 
statement (4) generalizes to: 


min 

— m + 1 j * « • ? Un , 
^n-m+l i ■ -^Otn) 


£ (^UfKiUi-Ffu^ + (iufKiUi-Fnui + RiaA 

1 = 1 V J 2=71 — 771+1 V / 


BjUi +■ Bj( u j + R jOLj) — 0 

i=l j=n— m+1 

R^Uj = 0 for i = n — m + 1, . . . , n 


(A2) 


The necessary conditions for the minimization problem given by statement (A2) are 


KjUj — Fj + Bf A = 0 


’V 


1 

CQ 

1 

fa 

i 

.Mi. 


0 


Rf(Fi-BfA) =0 


for i = 1, . . . , m — n 
for i — n~ m + 
for i = 


(A3) 


along with the first constraint equation given in statement (A2). Here A are the Lagrange multipliers corresponding 
to the first constraint equation in statement (A2), fa are the Lagrange multipliers for the last constraint equation 
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in statement ( 42 ), and the matrix K j for i — n - m + 1, . . . , n is the augmented stiffness matrix given by 


I 


K t 

R, 


, 

cT 


(AA) 


Using the second equation in necessary conditions (43), the condition R^K* = [0] for i = n — m + 1, . . . ,n, and 
the fact that the columns of are linearly independent, the last equation in necessary conditions (.43) implies that 
the Lagrange multipliers /Xf are equal to zero. The structural responses U* for i = 1, . . . , n — m, and (u ;,/!;) for 
i = n — m + 1, . . . , n are found by solving the first and second equations of necessary conditions (43), respectively, 
after the vector A is determined. 

The vector A must be determined simultaneously with the rigid body displacements a* for i = n — m + 
1, . . . ,n. This solution is accomplished by substituting the symbolic solutions for U j and u* into the compatibility 
constraint equation in statement (.42), and solving this equation simultaneously with the last equation in necessary 
conditions (43). The following matrix equation gives the resulting system: 



■ A ' 


' E?=r + Etn-m +1 [Bi I [0]]K-' [Ff I 0 T ] r - 


m+1 


~ B. n _ m+1 F n _ m+ i 

M 

&n—m - f 2 

= 

” B-^'_ m+2 F Il _ m+ 2 


Otn 


1 

e 

fa 

1 

i 


(45) 


where the matrix M is given by 


E”T 1 m B i K- 1 Bf+ j 

Er=„-m + i[B l |[0]]Kr 1 [B t |[0]] T 1 

1 

Rn-rn+l j 

1 

~ U n _ m +2Rn— m+2 j 

j “BnRn 

T 1 

“(Bn-m+lRn-m+l) j 

[0] 

[0] 

[0] 

T 1 

~ (Rn — m+2 Rn- m-f'i) j 

[0] j 

[0] 

[0] 

1 1 . 1 . 1 
1 : 1 : l *. I : 

! 1 II 

-(B„Rn) T | 

[0] 

[0] 

[0] . 


(46) 

Equations (45) and (46) are the generalizations of equations (9) and (10) for multiple substructures. 


Structural Optimization with Constraints Local to the Substructures 


As in the two-substructure example, only constraints that are local to the substructures (i.e. , that can be expressed 
using constraint functions of the form gi(vj,Ui) for i = 1 ,...,n) are considered in the structural optimization 
process. The steps utilized when there are multiple substructures are nearly identical to those performed for the 
two-substructure example. The simultaneous analysis and design formulation is a generalization of the problem of 
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statement (11) 


min 5Z^i(v t ) 

u n — m + 1 t 

Otn— m+ 1 r’l^ni 
Aj/in-m+lvipn) 

gi(vj,Uj)<0 for z = 1, . . . ,n — m 

gi(vj, + RiCti) < 0 for i — n — ra+l,...,n 

Ki(vj)Ut — Fj + B^A = 0 for i = l,...,n — m 


Ki(vt) 


Ut 


Fi-BfA’ 

.Mi. 


0 


= 0 for i = n- m+ l,...,n 


n— m 


yi BjUj + + ^ a j) _ o 

i=l j— n-m+l 

/it = 0 for i — n — m + l,...,n 
The Lagrangian function for this optimization problem is defined as 


(A7) 


n 


^ =51 Wt( v ») + 7 t T g«(vi, u i)+ ^rgiK-Ui + RiOi) 

t=l 1=71 — 771+1 


1=1 
n — m 


n — m n / 

+ ^ i/f (Ki(vj)Ut - F, + Bf A) + 5] if Kite) 

1=1 t=n— m +1 \ 


Mi 


Fi-BfA 

0 


(AS) 


i rt — lit ft i it 

+<57 (53 B i u i + 53 B j( u i +Rj«j) ) + 51 7, T Mt 

j=n— m +1 / i=n— m +1 


t=l 


and the following equations result from the optimality necessary conditions: 

dWi 


T^Si r dK,-Ui _ 0 r 

+7t a Vl +l/ ‘ 

, n ,T ^6* | T [ uf | o r ] T 

+~fi foT + Vi - 0 


dvi 


for i = 1 , . > . , n — m 
for i = n — m + l,...,n 


dv t i/vj 

7^ TJfF" + + 5 r Bj = 0 r for i = 






L [o] J 


n 




+ 6 T Hi = 0 T for i = n — m + 1, . . . , 
for i — n — m + 1, . . . , n 


ZWBT+ e -. t 

1=1 1=71 — 771+1 

[ 0 ] 


r Bp 

[ 0 ] 


- o' 


*7 K, 


l+i*9i 


+ T]f = 0 T for i = n — m + 


(*>) 


where p* is the number of degrees-of- freedom for the structural response of substructure i, and q{ is the number 
of rigid body modes in substructure i. This system of equations generalizes necessary conditions (13), and the 
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constraint Hi = 0 for i = n - m + 1, . . . ,n has been applied. The sensitivity derivatives of the structural response 
holding the interactions between the substructures fixed (i.e., fixed A) are defined as in definitons (15) 


dK,Vi 


for i = 1 , . . . , n — m 


^_ 1 gKj[uf |0 T ] T 

dui\ ~ i 


for i = n — m + 1 , . . . , n 


The Lagrange multipliers V{ are obtained from the equations in necessary conditions (A9) 


if | ^ 4- 6 T B^ 1 for i = 1, . . . , n - m 

~rf + S T Bj j nf K _1 for i = n - m + 1, . . . 
OUi 1 J 


(All) 


which are utilized with definitions (A10) to simplify the first two equations in necessary conditions (A9) 

dw i dgt , dgi dUj I \ 


dvi dV i dvi \J 

+5 r Bj = 0 T for i = l,...,n — m 


dWj T / dg, dgi duj 
d\i ^ I d\{ du i d\i 


{A12) 


i uv i \\ 


+6 T B x p + f p _ qT 

dvi A dvi j, 


= Cr for i — n — m + 1 , . . . , n 


As in the case of only two substructures, the final step is to determine values for 6 and r]i for i = n-m- hi, . . . , n. One 
equation necessary for determining these quantities is found by substituting expressions for from equations (All) 
into the sixth equation in necessary conditions (A9). The resulting equation becomes 

l ] T rE?.TBiK- 1 Bf + i:? = „- m+1 [B,!|0]]K- 1 [B,|[0]] T ' 

f?n-m+l m +iRn— m+l)^ _ 

~ _ /D T* \T (A13) 

_Vn—m+ lj L \&nK*n) J 

n—m o nr -i r "| T 

- £ - Y, I » T K . 1 B, ! [ 0 ] 

t=l 1 t=n— m-fl ^ J L J 

where the quantities i)i for i = n - m + 1, . . . ,n are defined by f)i = - [RfRi]~ X tj*. Combining equation (A13) 
with the fifth equation in necessary conditions (A9) yields the following equation that is solved to determine the 
values for 6 and 77 * for i = n — m + 1 , . . . , n 

* 1 [-Er,TB,K-'( 7 r|6t) T -Er„- m+1 [B,!io]]Kr 1 [ii r |g!o T ] T 

M Vn-m+l = (~fn-m+ 1 (A14) 
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Formulation of Decomposition as an Equilibrium Programming Problem 

The EP subproblems for n substructures are formulated by first defining approximate models for the displace- 
ments of the substructures that do not consider coupling between substructures. For example, linear approximate 
models may be used that have the form 


u " (v '> = 

"to)-*! 


(Vi - Vj) + Uj 


(Vi -Vi)+Uj 


for i = 1 , . . . , n — m 
for i = n — m + 1 , . . . , n 


(- 415 ) 


where the variables having caret accents are as defined previously, and the restricted derivatives are given by 
definitions (A10). These approximate models are utilized in formulating EP structural-sizing subproblems for the n 
substructures whose necessary conditions are the same as the equations given in reduced necessary conditions (A12). 
These structural-sizing subproblems are defined for the substructures having no rigid body modes (i.e., for 
i = 1 — m) as 


mmW 1 (v i )+6 T B l Vf(v t ) 

V l 


(-416) 


gi( v t,UiVi)) < 0 


The structural-sizing subproblems for the substructures that have rigid body modes (i.e., for i = n — m -f 1, . . . , n) 
are defined as 


min W 2 {vi) + 6 T Biuf{vi) + T]f ^ 

Vi 

gt(vi,u f(vi) + Rj«j) < 0 


Vi 


(-417) 


The Lagrange multipliers for in both statements (-416) and (-417) are represented by the vectors 7j which 
are needed to form the right side of equation (-414). The structural-sizing coordination subproblem given by 
equation (-414) determines the quantities 6 and fji for i = n — m 4- 1, . . . , n. The terms rfi are determined 
from rji — —HjH * 77^. The structural-response subproblems are given by the first two equations in necessary 
conditions (-43), and the structural- response coordination subproblem is given by equation (-45). 
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Figure 1 . Example wing structural finite element model 
decomposed into two substructures. 
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Abstract 


Quantum-mechanical optical model methods for calculating cross 
sections for the fragmentation of galactic cosmic ray nuclei by hy- 
drogen targets are presented. The fragmentation cross sections are 
calculated with an abrasion- ablation collision formalism . Elemental 
and isotopic cross sections are estimated and compared with mea- 
sured values for neon, sulfur, and calcium ions at incident energies 
between 400A MeV and 910 A MeV. Good agreement between theory 
and experiment is obtained. 

Introduction 

The fragmentation of galactic cosmic ray (GCR) nuclei in hydrogen targets is an important 
physical process in several areas of space radiation physics research. In astrophysics, it is crucial 
to understanding cosmic ray propagation and source abundances (ref. 1) because interstellar 
hydrogen is the major type of material encountered by GCR nuclei traveling through the 
universe. In studies of spacecraft shielding for interplanetary missions (ref. 2), hydrogen has 
been found to be the most effective GCR shield material per unit mass. In addition, hydrogen is 
a major constituent of human tissue. Therefore, accurate cross sections are needed for properly 
estimating GCR radiation exposures to critical body organs (ref. 3). 

Previously, cross-section predictions used in these studies have been obtained from semi- 
empirical formulations (refs. 4 to 7). The most commonly used formulation is the one by 
Silberberg and collaborators (ref. 5). The most accurate formulation appears to be a recent 
one by Webber and collaborators (ref. 6). None are based upon fundamental physics. All have 
numerous parameters that are adjusted as necessary to fit existing measurements. 

The production of fragments in peripheral, relativistic heavy ion collisions has been the 
subject of numerous theoretical and experimental investigations for about 2 decades. Many of 
these investigations were summarized in reviews published during this period (refs. 7 to 10). 
Early attempts to explain fragmentation used statistical models (refs. 11 and 12). These were 
followed by a two-step abrasion-ablation model (ref. 13), which was based upon earlier work by 
Serber in high-energy, inelastic nuclear collisions (ref. 14). 

The main shortcoming associated with the use of early abrasion-ablation models for nuclear 
fragmentation on hydrogen targets is the unrealistically large proton radius needed for the 
prefragment excitation energy estimate. This radius is dictated by the reliance on excess surface 
energy of the misshapen liquid drop as the only source of prefragment excitation. 

This shortcoming in the model can be rectified by considering an abrasion-ablation 
frictional-spectator-interaction (FSI) model where the abrasion stage is described by a quantum- 
mechanical optical model formalism and the ablation stage is modeled with cascade-evaporation 
techniques. There is no excess surface area energy. Instead, the prefragment excitation energy 
is assumed to be provided by FSI contributions from the abraded nucleons. This fragmentation 
model is proposed in this report. 

Abrasion- Ablation Models 

In an abrasion-ablation model, the projectile nuclei, moving at relativistic speeds, collide 
with stationary target nuclei. In the abrasion step (particle knockout), those portions of the 
nuclear volumes that overlap are sheared away by the collision. The remaining projectile piece, 
called a prefragment, continues its trajectory with essentially its precollision velocity. Because of 
the dynamics of the abrasion process, the prefragment is highly excited and subsequently decays 



by the emission of gamma radiation or nuclear particles. This step is the ablation stage. The 
resultant isotope is the nuclear fragment whose cross section is measured. The abrasion step is 
often formulated with methods obtained from quantum scattering theory (refs. 15 and 16) or 
with classical geometry arguments (refs. 13 and 17). The ablation step is typically modeled 
with compound nucleus decay (refs. 13 and 18) or combined cascade-evaporation (ref. 19) 
methods. Other approaches based upon nuclear Weiszacker-Williams methods (ref. 20) and 
nucleon-nucleon cascade plus statistical decay models (ref. 21) have also been proposed. 

Although abrasion-ablation fragmentation models have been quite successful in predicting 
fragment production cross sections, their predictive accuracy is hampered by the need to estimate 
the (unknown) prefragment excitation energy. Various models have been developed for this 
purpose (refs. 13, 15, 18, and 22). The most widely used excitation energy formalism (ref. 13) 
treats the fragmenting nucleus as a misshapen liquid drop whose excitation is given by the excess 
surface energy resulting from the abrasion step. Although this method worked fairly well for 
nucleus-nucleus fragmentations, its use in nucleus-hydrogen collisions, among other difficulties, 
required an artificially large proton radius (ref. 13). 

When it was recognized that additional excitation energy was required to improve the 
agreement between theory and experiment for nucleus-nucleus collisions, the concept of FSI 
energy was introduced (ref. 22). This concept is based upon the assumption that some abraded 
nucleons are scattered into rather than away from the prefragment, thereby depositing additional 
excitation energy. This concept significantly improved the agreement between theory and 
experiment. 

Over the past 10 years, we have formulated an optical model abrasion-ablation FSI descrip- 
tion of fragmentation in relativistic nucleus- nucleus collisions that is used to predict fragment 
production cross sections (refs. 16 and 23 to 42) and momentum distributions of the emitted 
fragments (refs. 43 through 47). In the present work, this fragmentation model is modified to 
make it applicable to nucleus-nucleon collisions. As previously discussed, the main shortcoming 
associated with the use of early abrasion-ablation models for nuclear fragmentation on hydrogen 
targets is the unrealistically large proton radius needed for the prefragment excitation energy 
estimate. This radius is dictated by the reliance on excess surface energy of the misshapen liquid 
drop as the only source of prefragment excitation. 

This shortcoming in the model can be rectified by considering the physics of the fragmentation 
process. For instance, a picture of overlapping nuclear volumes being sheared off may be 
reasonable for heavier nuclei colliding with each other, but it is not reasonable for a single 
nucleon striking another nucleus. Instead, a more reasonable physical picture involves individual 
collisions between the projectile constituents and the target proton. Some struck projectile 
nucleons exit the fragmenting nucleus without further interaction, and some interact one or more 
times with the remaining constituents before departing. The remaining nucleus (prefragment), in 
an excited state because of the energy deposited during the collision, then deexcites by particle- 
or gamma-emission processes. This picture is easily described by an abrasion-ablation-FSI model 
where the abrasion stage is described by a quantum-mechanical optical model formalism and 
the ablation stage is modeled with cascade-evaporation techniques. There is no excess surface 
area energy. Instead, the prefragment excitation energy is assumed to be provided by FSI 
contributions from the abraded nucleons. This fragmentation model is proposed in this report. 

Theory 

In the nucleus- nucleus optical potential formalism (ref. 29), the cross section for producing, 
by abrasion, a prefragment of charge Zpp and mass App is given by 

^abr (ZPF> j4pf) = J “ T(b)] n+ *[T(b)] ,4pF (1) 
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where 


T(b) = exp[— A t (r NN {e)I{h)] 


( 2 ) 


and 


/(b) = [27rJ5(e)] 3//2 J dz 0 J d 3 £ T p T (£ T ) J d A yp P (h + z 0 + y + £ r ) cxp[-y 2 /2B(e)} (3) 

The nuclear number densities pi{i ~ P or T) are obtained from the appropriate charge densities 
by an unfolding procedure (ref. 16). The constituent-averaged nucleon-nucleon cross sections 
a NN( e ) are given in reference 48. Values for the diffractive nucleon-nucleon scattering slope 
parameter B(e) are obtained from the parameterization in reference 49. 

In equation (1) a hypergeometric charge dispersion model is chosen to describe the distribu- 
tion of abraded nucleons. The model assumes that 2 out of Z projectile protons and n out of N 
projectile neutrons are abraded where 


N + Z = A P 


(4) 


v4pp = Ap — n — z 


(5) 


and 



denotes the usual binomial coefficient expression from probability theory. 


For nuclear collisions with hydrogen (proton) targets, the appropriate target number density 
to use is given by the Dirac delta function 


Pt(€t) -- 6(€t) ( 6 ) 

Inserting equation (6) into equation (3) yields 

I p (h) = [27 r5(e)] -3/ ' 2 J dz () J d 3 yp P {b + z 0 + y) exp[-y 2 /2B(e)} (7) 

With Ap = 1, equation (2) becomes 

T( b) = exp[— cr J vAr(e)/p(b)] (8) 

The nucleus- hydrogen abrasion cross sections are calculated with equations (1), (7), and (8). 
Prefragment excitation energies are estimated from the FSI energy contribution 


£exc = Epsi 


(9) 


which is calculated with the model of Rasmussen (ref. 22). With this model, the rate of energy 
transfer to the prefragment is 


dE _ E_ 
dx 4A 


(10) 


where 


P a NN 


(°NN ~ 



( 11 ) 
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yields 

— = -12.75 McV/fm (12) 

dx 

If a spherical nucleus of uniform density is assumed, the average energy deposited per 
interaction is 

{E F si) « 10.2A 1 / 3 MeV (13) 

Therefore, the abrasion cross section for a prefragment species (Zpp,App) which has 
undergone q frictional spectator interactions is 

a. A \ )T (Zpp, TpF,<?) = _ ^esc) 9 (^esc) n+2 ~ 9cr abr(^PF)^PF) ( 14 ) 

where 0 < q < n + z, and P esc is the probability that an abraded nucleon escapes without 
undergoing any frictional spectator interactions (ref. 34). In this report, the choice of Pesc = 0-5 
follows from the original work of Rasmussen (ref. 22). Such a value assumes that the nuclear 
surface has no curvature, and this value should be reasonably correct for heavy nuclei. For 
lighter nuclei, the surface can exhibit significant curvature such that the value of P cs c can be 
larger than 0.5. Methods for estimating P esc when nuclear surface curvature is considered have 
been formulated by Vary and collaborators (ref. 50). 

Depending upon the magnitude of its excitation energy, the prefragment will decay by 
emitting nucleons, composites, and gamma rays. The probability a,j(g ) that a prefragment 
species j, which has undergone q frictional spectator interactions, deexcites to produce a 
particular final fragment of type i is obtained with the EVA-3 Monte Carlo cascade-evaporation 
computer code (ref. 19). Therefore, the final hadronic cross section for production of the type i 
isotope is obtained from 


n+z 

0nuc(^7 •> Aj) ~ EE OLij{q)(T & bAZj,Aj,q) (15) 

j <7=0 

where the summation over j accounts for contributions from different prefragment isotopes j, 
and the summation over q accounts for the effects of different FSI excitation energies. Finally, 
the elemental production cross sections are obtained by summing all isotopes of a given element 
according to 

^nuc {%i) = ^ ^ fl'nuc ( Zj ? Aj ) (16) 

A,: 


Results 

Figures 1 and 2 show isotope production cross sections obtained with equation (15) for ’* 2 S 
beams at 400A MeV fragmenting in hydrogen targets. The figures also show recently reported 
experimental results (ref. 51). For clarity, the experimental error bars are not plotted. The 32 S 
nuclear density used in the calculation was a Woods-Saxon form with skin thickness and half- 
density radius obtained from reference 48. The agreement between theory and experiment is 
quite good, especially considering that no arbitrary parameters are in the theory. Quantitatively, 
a distribution analysis of the cross-section differences between theory and experiment finds that 
32 percent agree within the experimental uncertainties, 50 percent agree within a 25-percent 
difference, nearly 75 percent agree within a 50-percent difference, and over 82 percent agree 
within a factor of 2. 

Elemental production cross-section predictions obtained from equation (16) are displayed 
in figures 3 to 8 for 2,) Ne beams at 400A MeV and 910A MeV and for 32 S and 40 Ca beams at 
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400A MeV and 800A MeV incident kinetic energies colliding with hydrogen targets. The nuclear 
densities used in the calculations were Woods-Saxon forms with skin thicknesses and half-density 
radii again obtained from reference 48. These experimental data were taken from reference 51. 
Overall, the agreement between theory and experiment is good, although the theory tends to 
predict values that are slightly larger than the reported measurements. 

• Theory 
□ Experimental F 
2 ■ Experimental Na 

jq f ♦ Experimental Ai 
E o Experimental P -/* 


12 14 16 18 20 22 24 26 28 30 32 
Fragment mass 

Figure 1. Isotope production cross sections for 400 A MeV 32 S fragmentation in hydrogen targets for isotopes of P, 
AL, Na, and F fragments. 

• Theory 
o Experimental O 
-I q2 a Experimental Ne 
_o _ o Experimental Mg ■r'A 

E □ Experimental Si [ \ D 

!:: a/Vx 


12 14 16 18 20 22 24 26 28 30 32 
Fragment mass 


Figure 2. Isotope production cross sections for 400A MeV 32 S fragmentation in hydrogen targets for isotopes of Si, 
Mg, Ne, and O fragments. 



Fragment charge 


Figure 3. Element production cross sections for 400A MeV 2{) Ne fragmentation in hydrogen targets. 
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Fragment charge 

Figure 4. Element production cross sections for 910A MeV M Ne fragmentation in hydrogen targets. 



4 8 12 16 

Fragment charge 


Figure 5. Element production cross sections for 400A MeV 32 S fragmentation in hydrogen targets. 



4 8 12 16 

Fragment charge 


Figure 6. Element production cross sections for 800A MeV 32 S fragmentation in hydrogen targets. 
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Fragment charge 

Figure 7. Element production cross sections for 400A MeV 40 Ca fragmentation in hydrogen targets. 



Fragment charge 

Figure 8. Element production cross sections for 800A MeV ll 'Ca fragmentation in hydrogen targets. 

Concluding Remarks 

A simple, yet accurate, optical potential abrasion-ablation fragmentation model has been 
developed for use in studies of galactic cosmic ray breakup on hydrogen targets. The model 
has no arbitrarily adjusted parameters. Model predictions have good agreement with recent 
laboratory measurements of elemental and isotopic production cross sections for the fragmenting 
of neon, sulfur, and calcium beams on hydrogen targets. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
October 28, 1993 
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